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ABSTRACT 


i(jk  semi-atomistic  model  is  proposed  for  the  non-ductile 
failure  of  materials  based  upon  an  analysis  of  static  and  dynamic 
stress  concentrations  in  uniaxial  fiber-matrix  composite  materials. 
In  this  model  the  finite  spacing  of  bonds  is  accounted  for  in  the 
direction  of  a  crack  or  flaw,  while  the  material  is  considered  as 
continuous  in  the  perpendicular  direction.  An  elementary,  two- 
dimensional  analysis  based  on  linear  force-displacement  laws  shows, 
that  the  stress  concentration  is  dependent  only  upon  the  digital 
quantity  of  the  number  of  adjacent,  broken  filaments  or  bonds,  and 
is  independent  of  the  spacing  and  of  the  relative  elastic  properties 
of  the  fiber  and  matrix.^  The  form  of  the  static  stress  concentra¬ 
tion  factor 


where  n  is  the  number  of  adjacent  broken  bonds,  c  is  the  crack  depth 
and  d  is  the  bond  spacing,  is  formally  equivalent  to  that  obtained 
using  relations  given  by  Griffith  (Ref.  1)  and  Charles  (Ref.  3) . 

The  results  of  the  elementary  analysis  are  of  sufficient 
significance  to  warrant  extension  of  the  theory  to  non-linear  force 
elongation,  and  to  account  for  temperature,  stress  corrosion  and 
static  fatigue  effects.  In  addition,  experimental  verification  of 
this  theory  should  be  undertaken  using  both  covalent  and  ionically 
bonded  materials ,  as  well  as  metals  under  conditions  where  non 
ductile  failure  occurs. 
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I.  INTRODUCTION 


The  technical  strength  of  non-ductile  materials  has  become  of 
increasing  interest  because  of  their  use  in  fibers  (glass,  elemental 
boron,  oxides,  carbides,  etc.)  and  thermosetting  polymers  (epoxies, 
silicones,  polybenzimidazoles,  etc.)  that  form  the  constituents  of 
many  high  performance  composites  used  in  advanced  structural  applica¬ 
tions.  Further,  many  methods  of  improving  the  tensile  strength  of 
metallic  materials  do  so  by  reducing  dislocation  mobility  and  tend 
towards  generating  non-ductile  fracture  modes.  Also,  the  use  of 
refractory  oxide  ceramics ,  carbides ,  graphite  and  carbon  in  bulk 
shapes  for  structural  applications  forces  attention  upon  their 
essentially  non-ductile  failure  phenomena. 

Other  reasons  for  interest  in  the  strength  and  failure  of  non- 
ductile  materials  are  found  in  cryogenic  applications.  In  these 
applications,  elastomeric  or  thermoplastic  materials  are  often 
forced  to  operate  at  temperatures  below  the  "glass  transition", 
where  they  exhibit  brittle,  non-ductile  fracture  associated  with 
catastrophic  crack  propagation. 


II.  INTRINSIC  AND  TECHNICAL  STRENGTH  OF  MATERIALS 


Solids  derive  their  resistance  to  deformation  from  three 
types  of  interatomic  bonds: 
metallic, 
ionic,  and 
covalent. 

A  schematic  representation  of  the  force-deformation  characteristic 
inherent  in  the  interatomic  bond  is  shown  in  Figure  1. 

Conventional  materials  of  construction  are  predominantly 
metallic.  Metals  exhibit  dislocation  mobility,  which  allows  the 
relief  of  stress  concentration  by  what  is,  in  effect,  a  premature 
but  non-catastrophic  (i.e.,  "self-healing")  failure  mode.  A  general 
term  for  this  type  of  materials  response  to  stress  is  "ductility". 
Even  metallic  materials,  however,  may  exhibit  non-ductile  fracture 
where  impurities,  work  hardening  and  cyclic  fatigue  (dislocation 
pile-up) ,  or  peculiar  crystallographic  properties  (beryllium) 
restrict  dislocation  mobility. 

Other  types  of  materials  derive  their  cohesion  from  electro¬ 
static  forces  (ionic  bonds)  and  from  quantum  mechanical  electron 
exchange  forces  (covalent  bonds) .  The  latter  bonds  are  typical 
for  ceramics,  glasses,  etc.  Because  covalent  bonding  forces  are 
space  oriented  they  usually  prevent  dislocation  motion,  and  therefore, 
in  principle,  can  provide  materials  of  much  higher  strength. 
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However,  because  of  the  absence  of  the  "ductile"  stress  relief 
mechanism  possible  for  metallic  (and  some  ionic)  bonds,  covalently 
bonded  materials  are  subject  to  other  premature  failure  modes  that 
are  associated  with  either  static  'or  dynamic  stress  concentrations 
(brittle  fracture) . 

The  observed  tensile  strength  of  non-ductile  materials  is 
orders  of  magnitude  lower  than  the  intrinsic  interatomic  cohesion 
that  can  be  derived,  for  instance,  from  their  heat  of  sublimation. 
Following  original  work  of  Griffith,  Inglis,  and  others  (Refs.  1-3) 
explanations  for  this  deficiency  have  been  given  in  terms  of  ob¬ 
served  or  postulated  flaws  that  are  responsible  for  microscopic 
stress  concentrations  in  a  seemingly  uniformly  stressed  body. 

One  unsatisfactory  aspect  of  the  classical  flaw  theory  is  its 
dependence  upon  an  isotropic,  continuous  medium  for  calculation  of 
the  stress  at  the  crack  tip  and  the  necessity  to  postulate  a  finite 
crack  tip  radius  to  obtain  finite  values  for  stress  concentration 
factors.  Equally  unsatisfactory  are  treatments  involving  a  free- 
surface  energy  or  "surface  tension"  and  related  potential  energy 
postulates  to  determine  critical  flaw  size.  The  first  approach 
obviously  violates  the  discrete  nature  of  molecular  arrangements 
that  must  be  considered  in  an  atomistic  view  of  the  crack  tip 
geometry.  The  second  approach  deals  with  the  static  situation 
prevailing  at  some  time  after  the  failure  has  occurred,  and  is 
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incapable  of  accounting  for  the  detailed  dynamic  mechanisms  attend¬ 
ing  the  straining  and  failure  of  interatomic  bonds.  (A  portion  of 
the  failure  energy  must  leave  the  failure  domain  in  the  form  of 
thermoelastic  wavelets  or  phonons.)' 

An  attempt  to  account  for  finite  bond  spacing  has  been  presented 
by  Elliot  (Ref.  4) .  Here  the  material  was  considered  as  a  Hookean 
elastic  continuum  with  the  exception  of  two  monatomic  layers  immedi¬ 
ately  adjoining  the  crack  area.  This  concept  allows  accounting  for 
the  non-linear  force-elongation  characteristics  of  highly  strained 
bonds,  but  suffers  from  the  transition  to  continuum  mechanics  at  an 
atomic  distance  removed  from  the  crack.  For  instance,  analysis  of 
cohesive  force  distribution  based  on  Elliot's  model  yields  large 
stress  gradients  in  distances  smaller  than  the  interatomic  spacing. 
This  result  is  difficult  to  reconcile  with  the  discrete  nature  of 
materials . 

A  new  approach,  suggested  by  the  results  of  a  recent  study  on 
crack  propagation  in  filamentary  materials  with  partially  broken 
filaments  (Ref.  5) ,  avoids  the  necessity  of  postulating  crack  tip 
radii,  but  retains  an  essentially  continuum  mechanical  approach 
suitable  for  analytical  treatment.  As  applied  to  interatomic 
dimensions  this  approach  provides  a  semi-atomistic  model  in  the 
sense  that  finite  spacing  of  bonds  is  accounted  for  in  the  direction 
of  the  crack,  while  the  material  is  considered  as  continuous  in  the 
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direction  perpendicular  to  the  crack.  This  concept  is  developed 
here  for  the  prediction  of  flaw  effects  on  technical  strength ,  and 
is  correlated  with  published  experimental  data*  Extension  of  this 
theory  to  non-linear  force-elongation,  and  to  account  for  tempera¬ 
ture  effects,  stress  corrosion,  static  fatigue,  etc.,  is  possible 
but  has  not  been  carried  out  in  detail* 


in  principle, 


III. 


FRACTURE  ANALYSIS 


Hedgepeth  (Ref.  5)  has  treated  the  case  of  a  one-dimensional 
array  of  parallel  elastic  filaments  in  a  matrix  with  elastic  shear¬ 
ing  stiffness,  and  with  a  crack  extending  over  n  filaments,  as 
shown  in  Figure  2.  The  array  is  assumed  to  be  subject  to  uniaxial 
tension  in  a  direction  normal  to  the  crack.  The  static  stress  con¬ 
centration  factor,  Kn,  for  the  first  unbroken  filament  is  given  by 
the  recursion  formula 

K  -  —  -  6  •  3  (2n  +  2) 
n  ~  3  .  5  •  7  (2n  +  1)  1 


Figure  3  shows  the  relation  between  the  number  of  broken  filaments, 
n,  and  the  stress  concentration  factor,  as  given  in  Equation  1. 

Kn  diverges  for  large  n,  but  can  be  represented  with  high  accuracy 
for  n  >  10  by 


Since  exp  converges  rapidly  to  unity  for  large  n.  Equation  2 

can  also  be  approximated  by 


*  The  derivation  of  this  expression  is  due  to  G.  Schindler  and 
is  given  in  Appendix  A. 
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It  is  significant  that,  within  the  limitations  of  the  analysis, 
the  stress  concentration  is  independent  of  the  spacing  and  of  the 
elastic  properties  of  fiber  and  matrix,  and  is  dependent  only  on 
a  digital  quantity,  i.e.,  the  number  of  adjacent  broken  filaments. 

It  is  further  interesting  that  the  form  of  the  stress  con¬ 
centration  factor  given  in  Equation  2a  is  formally  equivalent  to 
that  which  can  be  obtained  by  using  relations  given  by  Griffith 
(Ref.  1)  and  by  Charles  (Ref.  3)  for  the  critical  stress  at  failure 
°cr'  and  for  the  intrinsic  strength,  F,  respectively: 


where  S  is  the  surface  energy,  E  is  the  Young's  modulus,  c  is 
the  crack  depth,  and  a  is  the  interatomic  spacing. 

Combining  Equations  (3)  and  (4)  yields 


In  addition  to  the  static  stress  concentration  described  by 
Equation  1,  dynamic  overload  stresses  occur  due  to  the  elastic 
waves  generated  by  the  sudden  fracture  of  one  or  several  filaments. 
The  dynamic  overload,  n^,  due  to  a  sudden  fracture  of  k  filaments 


is  given  by  Hedgepeth  (Ref.  5)  for  1,  2/  3  filaments  as  1.15, 

1.19  and  1.20,  respectively,  with  a  limiting  value  of  =  1.27 

for  the  simultaneous  failure  of  many  filaments. 

This  analysis  assumes  that  k  filaments  fail  simultaneously, 
but  that  prior  to  the  failure  all  filaments  were  intact.  The 
results  of  this  analysis  can  also  be  used  as  a  first  approximation 
to  the  case  where  pre-existing  failures  have  generated  a  static 
crack,  and  where  the  instantaneous  failure  of  k  additional  fila¬ 
ments  produce  a  dynamic  stress  concentration. 

Consider  now  a  non-ductile  material  with  a  "flaw"  in  the  form 
of  a  plane  region  perpendicular  to  the  direction  of  applied  stress, 
as  shown  in  Figure  4.  The  flaw  may  either  be  at  the  surface  or  in 
the  interior  of  the  body. 

We  may  now  use  the  concept  of  "filaments"  for  the  arrays  of 
molecular  bonds  in  the  direction  of  the  applied  stress.  A  plane 
section  parallel  to  the  axis  of  stress,  and  passing  through  the 

diameter  of  the  flaw,  may  then  exhibit  a  picture  very  much 
like  the  two-dimensional  array  of  filaments  shown  in  Figure  2. 

The  stress  concentration  prevailing  in  the  bonds  adjacent  to  the 
flaw  edge  is  given  by  Equation  1  or  2 ,  and  is  seen  to  depend  solely 
on  the  number  of  broken  bonds,  n  =  2c/d.  In  this  equation  2c  is 
the  maximum  diameter  of  the  interior  flaw  (or,  c  is  the  depth  of  a 
surface  flaw),  and  d  is  the  spacing  of  the  broken  bonds,  measured 
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in  the  plane  of  the  flaw.  With  this  semi-atomistic  model,  the 
critical  stress  level  at  which  crack  propagation  will  occur  becomes, 
for  large  n. 


(6) 


It  is  seen  that,  according  to  this  analysis,  failure  will 
occur  at  a  fraction  of  the  intrinsic  strength,  and  that  this 
fraction  is  dependent  both  on  the  number  of  broken  bonds  in  the 
pre-existing  flaw,  and  on  the  nature  of  the  micro-dynamic  stress 
concentration  caused  by  sudden  additional  bond  fracture. 
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IV.  DISCUSSION  OF  MICRO-DYNAMIC  EFFECTS 


We  may  distinguish  three  mechanisms  of  crack-propagation, 
depending  upon  severity  of  the  dynamic  effect: 

a)  Quasi-static  (Creeping)  Crack  Propagation 

Crack  propagation  is  thermochemically  activated,  i.e., 
the  elastic  wave  (phonon-packet)  emitted  by  the  failure 
of  a  single  bond  is  insufficient  to  cause  immediate 
failure  of  additional  bonds.  The  static  stress  concen¬ 
tration  is  high  enough,  however,  that  after  sufficient 
elapsed  time,  failure  will  occur  by  chemical  attack  of 
adsorbed  atmospheric  constituents  (particularly  of  water 
in  the  case  of  silica  glasses)  and,  possibly,  by  random 
thermal  motion.  This  process  is  similar  to  that  which 
can  be  postulated  for  surface  solubility  and  sublimation 
of  solids  except  that  the  bond  failure  probability  here 
is  enhanced  by  the  mechanically  induced  bond  strain  at 
the  crack  tip.  These  effects  will  reduce  the  intrinsic 
bond  strength,  but,  in  this  case,  the  "dynamic  stress 
concentration"  can  be  neglected, 
b)  Transition 

There  is  a  threshold  case  where  the  transient  stress 
concentration  of  the  bond  failure  process  is  sufficient 
to  cause  failure  of  the  adjacent  bond  before  the  dynamic 
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stress  has  dissipated ,  but  with  sufficient  time  delay 
to  prevent  the  catastrophic  interaction  of  phonons 
emitted  by  consecutive  failures •  The  speed  of  crack 
propagation  will  be  expected  to  be  small  compared  to 
the  velocity  of  sound  in  the  material.  Dynamic  over- 
load  factors  of  the  order  of  those  given  in  Reference  5 
for  single  instantaneous  bond  failure  will  be  expected. 

)  Catastrophic  Crack  Propagation 

Crack  propagation  is  "instantaneous"  in  the  sense  that 
the  dynamic  effect  of  a  large  number  of  bond  failures 
will  accumulate  causing  the  bond  at  the  crack  tip  to 
fail  without  thermal  delay.  Thus,  a  shock  wave  will  be 
generated  that  will  be  expected  to  propagate  at  sonic 
velocity. 


-V.  EXPERIMENTAL  CORRELATION 


Experiments  concerning  the  tensile  strength  of  glass  with 
flaws  or  initial  cracks  of  various  sizes  and  origins  have  been 
reported  by  many  investigators .  The  absolute  values  of  measured 
strength  vary  over  a  wide  range  depending  on  conditions  of  test, 
size  and  shape  of  test  specimens,  atmosphere,  aging,  previous  heat 
treatments  of  glass  (temper) ,  etc.  It  is  found  by  many  investigators, 
however,  that  the  measured  strength  is  inversely  proportional  to  the 
square  root  of  the  flaw  size,  as  predicted  by  Equation  6.  Also,  a 
marked  dependency  of  failure  strength  with  rate  of  load  application, 
"static  fatigue",  and  slow  growth  of  cracks  that  are  originally 
smaller  than  the  critical  size,  has  been  well  established.  These 
observations  correlate  qualitatively  with  the  dynamic  failure  model 
discussed  here. 

The  theoretical  critical  bond  failure  number  n,  given  by 
Equation  6,  is  plotted  vs  a  /F  in  Figure  5,  using  Hedgepeth's 

V/  JL 

"two-dimensional"  dynamic  factors  of  1.15  and  1.27  for  transition 
and  catastrophic  crack  propagation,  respectively,  and  using  1.0  for 
the  quasi-static  case  of  "creeping"  crack  propagation.  Also  plotted 
in  Figure  5  are  data  from  Reference  1  (for  relatively  large  cracks) 
and  Reference  4  (for  small  surface  abrasions)  for  comparison  between 
the  present  theory  and  experiments. 
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For  the  purpose  of  this  comparison,  several  assumptions  were 


made : 

1.  The  bond  spacing  d  was  assumed  to  be  the  inverse  square 
root  of  the  number  of  Si-0  bonds/unit  area  in  fused 
silica  (Ref.  3)  yielding 

d  =  5.4  A  =  2.1  x  10-8  in. 

We  note  that  this  value  is  considerably  larger  than  the 
Si-0  bond  distance 

O  —ft 

a  .  =  1.6  A  =  .63  ^  10  in. 

Si  0 

The  reason  for  this  difference  becomes  evident  if  the 
open  network-structure  of  silica  shown  schematically  in 
Figure  6  is  considered. 

2.  The  intrinsic  strength  was  assumed  to  be  2 . 5  x  10^  psi 

for  glass  tested  in  normal  atmosphere  (tests  by  Griffith, 

Ref.  1)  1.8  x  106  psi  and  3.2  x  106  psi  for  tests  in  water 

and  in  liquid  nitrogen,  respectively  (tests  by  Mould  and 

£ 

Southwick,  Ref.  6).  The  value  of  2.5  x  10  psi  correlates 
well  with  experimental  values  obtained  from  thin,  presumably 
flaw-free  filaments  (Ref.  3) ,  and  with  the  theoretical 
strength  obtained  from  Equation  4,  if  the  bond  distance,  a, 
is  replaced  by  the  bond  spacing  d  =  5.4  A.  The  increase 
of  intrinsic  strength  due  to  immersion  in  liquid  nitrogen 
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may  be  correlated  to  two  factors :  The  lower  temperature 
reduces  thermal  activation  of  the  failure  mechanism  itself. 
Also,  the  lower  temperature  will  reduce  chemical  attack  of 
the  fracture  tip  with  attendant  lowering  of  the  residual 
intrinsic  bond  strength  such  as  observed  in  the  tests 
conducted  immersed  in  water. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  semi-atomistic  model  of  non-ductile  failure  appears  to 
predict  the  essential  features  of  failure  observed  in  glassy  ma¬ 
terials.  The  analysis  here  has  been  restricted  to  a  two-dimensional 
model  and  has  not  accounted  for  the  detailed  mechanisms  attending 
an  intermolecular  bond  failure.  This  refinement  is,  however, 
possible  in  principle  and  should  yield  a  more  complete  understanding 
of  failure  mechanics  on  an  atomic  level.  Two  particular  subjects 
should  be  considered: 

-  A  theoretical  treatment  of  the  dynamic  failure  phenomena, 
including  consideration  of  thermal  vibrations,  and 

-  An  expansion  of  Hedgepeth's  work  to  include  three- 
dimensional  arrays  of  bonds,  as  well  as  network-structures 
of  the  type  shown  in  Figure  5,  and  to  include  realistic 
non-linear  force  displacement  characteristics  of  the 
intermolecular  bond. 

Well  controlled  static  and  dynamic  fracture  experiments  with 
non-ductile  materials  of  current  technical  interest,  including 
artificially  induced  flaws  in  the  test  specimens  should  be  conducted 
to  correlate  refined  theoretical  predictions.  Such  experiments 
should  be  conducted  first  in  inert  atmospheres  to  eliminate  the 
obscuring  effects  of  chemical  crack  tip  attack.  Later,  both 
theoretical  and  experimental  work  should  be  devoted  to  a  quantitative 
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understanding  of  "stress-corrosion"  as  it  affects  catastrophic 
and  creeping  crack  growth. 

Candidate  materials  for  such  experimental  studies  are  refractory 
metals,  thermosetting  polymers,  thermoplastics  and  elastomers  at 
temperatures  below  the  glass  transition,  and  "exotic"  polycrystalline 
or  amorphous  materials  such  as  boron,  beryllium,  carbon  and  graphites 
in  either  bulk  or  fiber  form. 
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FIG.  h  INTRINSIC  STRENGTH  OF  MATERIALS 
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FIG.  6 !  SCHEMATIC  REPRESENTATION  OF  GLASSY 
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APPENDIX 


Reference  5  determined  the  stress  concentration  factor,  KR/ 
for  a  plane  filamentary  structure  in  the  static  equilibrium  in 
the  case  of  n  broken  filaments  as  follows: 


4  •  6  •  8 
3-5-7 


( 2n  +  2) 
(2n  +  1) 


( A-l) 


For  large  n,  expression  A-l  is  inconvenient  to  evaluate 
Equation  A-l  can  be  written 


K  =  2 

n 


2n  n!  (n+1)  ! 
( 2n+l) ! 


According  to  Reference  7,  the  quantity  n!  can  be  expressed  in 


the  form 


n!  =  v  2  it  n 


-n+A. 


(A-3) 


where 


-  f  ± 

J  x 


sin  2irrx  <Jx 


x  2 
r=l 


( A-4 ) 


Inserting  (A-3)  into  (A-2) ,  one  obtains 


K  =  ±  v  irn 

n  2 


(l+i 

.  V _ n 


(  An+An+l"A2n+l )  (A-5 


1  \  2n+: 


Forming  the  natural  logarithm  of  this  expression,  it  follows  that 


A-l 


109  K„  -  i  109(f)  +  (n+|)  log  (1+  I)-  (2n+|)  log  (l+^) 


(An  +  An+1  '  A 


2n+l 


(A-6) 


If  n  i  1,  the  expressions  log  (1+— )  an^  1 og  (l+— j  can  be  written 


m  a  power  series 


log  (  1+i  )  =  i 
v  n  /  n 


1-  +  I 


1  - 
3  +  ’ 


=  I  (-l)r_1  I  .  /  1  \  r 


(I) 


1o9(1+^)  = 


.1-1.1  .1  1 _  T 

2n  2  (2n)  2  3  (2n)  3 


2  (-1) 
r=l 


After  some  transformations  we  obtain 


ni) 109  (1+~) 


2n+|j  log  ^l+_i  j 


1_  +  I 


1 _  .  B 

(2n)2  '  n 


Inserting  (A-7)  into  (A-6)  yields 


I09  Kn  -  I  109(f)  *  f 


+  I 


(1ST2  Bn  +  (A"  +  An+1  '  A2n+1 


(A-8) 


According  to  Equation  (A-7) ,  the  series  for  Bn  can  be  majorized 


as  follows 


A-2 


B 

4 

oo 

z  1  • 

[4  (i)r  - 

r 

( S-) 

64  \(b) 

n 

r=0 

\  n  / 

\  2n  / 

r=0  \ n  / 

1  - 


n 


Equation  (A-4)  yields  after  integration  by  parts: 


_1  £ 


A  =  -  0 

n  n  r=l  4ii2  .  r2 


-  4 


i 


n 


Now,  according  to  Reference  7, 


1 

Z  ” 
r4 

r=l 


1  i  sin2Trrx 


x3  r=]_  ( 27rr)  3 


dx 


(A 


Hence : 


0  ^  Z 

r=l 


^  Z  “7 
3  2 
r  r=l  r 


A 


n 


1_ 

12 


.  -  -  C 


n 


1-  .  i  -  4 

12  n 


i  1 


sin2iTrx 


(A- 


dx 


r=l  (27rr) 
n 

The  sum  (k^  +  An+1  -  A2n+l)  Can  be  exPressed  aS  follows: 
An  +  An+1  -  A2n+1  : 

^  Cn  +  Cn+i  “  C2n+1  ^ 


1_ 

12 


I  +  Jk _ 

n  n+1  2n+l 


(A-9 ) 


-10) 


ID 


(A- 12 ) 


(A-13) 


A-5 


The  remainder  can  be  estimated  by 


Cn  +  Cn+1  "  C2n+1  ^  4  •  3  J 


-  £  - ^  q  v  ^ 

x3  r,1  (2»r)3 


r=l  r 


1  4 


dx  = 


(A-14 ) 


It  is  further  for  n  >  0: 


(4) 


(A-15) 


_1 _ 

2n+l 


Hence,  Equation  (A-13)  reads 


A  i  A  _  A 

n  n+1  2n+l 


1  1  _  1_ 
8  '  n  12 


(>4) 


1  1 
4  '  7" 


1  +  ~  n' 
2n/ J 


+  Si+1  *“2n+l 


(A-16 ) 


Equation  (A-8)  reads 


or 

log  K  =  -  log  /  — )  +  -  •  -  +  A  (A-17) 

n  2  ^  \  4  j  8  n 


The  remainder  A  can  be  estimated: 


Equations  (A-9)  and  A-14)  yield 


A  .  1  1 

i 

4 

u- 

5 

+  I. 

A  -  . 

n 

8 

(4) 

T 

48 

8  TT 

For  n  i  10  we  have 


1  A  ^  1 

1 

+  _s 

1_  j 

S' 

1 

1  A  100 

2(1 -  h.) 
\  10 ; 

1  48 

8  TT  j 

100 

Hence,  for  n  i  10,  log  Kn  can  be  expressed  with  sufficient 
accuracy  by 


A-5 


log  K. 


2  log 


(f) 


*1';  =  °-01 


(A-20) 


Since 


-0.01  +  0.01 

0  ^  e  e  =  1.01005  . 


1.011 


we  obtain  the  approximation  formula 


n  2 


_  5.1 

un  e  8  n 


(A-21) 


which  gives  the  exact  value  for  Kn  within  an  accuracy  of  less 


than  1.1%,  if  n  2  10.  It  is  clear  that  tends  to  infinity  like 
the  square  root  of  n. 


A-6 
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